Power-law behavior observed in j^^-distributions and its 
implications in relativistic heavy-ion collisions 
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Preconception-free analyses of the inclusive invariant transverse-momentum dis- 
tribution data taken from the measurements of Au-|-Au collisions at ^sjvAf = 130 



, GeV and ^/snn = 200 GeV have been performed. It is observed that the distribu- 

I tions exhibit for ^ 2 GeV/c remarkably good power-law behavior (p^-scaling) 

O ' 

^SJ ■ with general regularities. This power-law behavior leads us in particular to recognize 

^ I 

. that the concept of centrality, albeit its simple appearance, is rather complex; its un- 

P 

derlying geometrical structure has to be understood in terms of fractal dimensions. 
Experimental evidences and theoretical arguments are given which show that the 

^ : 

2 ■ observed striking features are mainly due to geometry and self-organized criticality. 

A simple model is proposed which approximately reproduces the above-mentioned 
^ I data for the "suppression" without any adjustable parameter. Further heavy-ion 

^ collision experiments are suggested. 

PACS numbers: 25.75.-q, 25.75.Dw, 05.65. +b, 13.85.Hd, 13.87.Ce 

Inclusive invariant pT-distributions for charged hadrons in Au+Au collisions at a/s/vw" = 
130 and 200 GeV have been measured and published by STAR and by PHENIX 3] over 
a broad range of centrality. Such »T-distributions for neutral pions are also given by PHENIX 
0] . All these experiments Q, Q, 0,1^] show that the hadron yields differ appreciably at high 
and medium px in central collisions relative to peripheral collisions and to the nucleon- 
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FIG. 1: Inclusive invariant pT-distribution data 
white indicate 130 and 200 GeV respectively. 

nucleon reference. What do these observations, usually known as "suppression" P], tell 
us? Are they related to the yet-to-be-found "quark-gluon-plasma (QGP)"? If yes, how? In 
order to obtain an unbiased physical picture to start with, we begin with preconception-free 
data-analyses. We then summarize the results, and discuss their implications. 

In the first part of this paper, we report on the result of such analyses. Within the 
measured kinematical region < pt < 12 GeV/c, ^/sNN = 130 GeV and ^/sNN = 200 GeV, 
the inclusive invariant pr-distributions of {h^ + h^)/2 for centrality-selected Au+Au, and 
those for p+p interactions exhibit power-law behavior for Pt > 2 GeV/c (cf. Fig. 1). The 
results can be summarized as follows. 
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where the power-indices (the A's) are positive real numbers, and Pc characterizes the 
centrality-bins (in percentile) which stand for the different degrees of departure from the 
most central collision. The experimental values of the A's obtained from the STAR data 



at ^/SNN = 130 GeV and 200 GeV for different pc-bins are shown in Fig. 2. The 



results from the PHENIX data 



^ (will be reported in a more extensive paper elsewhere 
show similar characteristic properties. 
As we can see in Fig. 2, for a given Pc-hin, the A- value at ^/sNN = 200 GeV coincides 
with the corresponding one at a/J/vtv = 130 GeV. Furthermore, the A-value for the most 
peripheral {pc — > 100%) case is very much the same as App's. Note that the A-values increase 
from the most peripheral value, XauAu{Pc 100%) ~ App, with decreasing pc to the A-value 
for the most central (center-on-certer) collision {p^ 0%) in a monotonous manner. In this 
connection it is useful to consider the ratio of both sides of Eqs. (1) and (2): 

(fN/pTdpTdr][^=o{pT\Au + Au; pc) -A(p,) 



o^Pt'"''- (3) 



d"^ N / pTdpTdri\r,=o{pT\p + P', inel.orNSD) 
The quotient, which is a completely experimental quantity, on the left-hand-side of this 
equation, will hereafter be referred to as Q{pt,Pc)] and the values of the exponent on the 
right-hand-side, 

HPc) = >^AuAu{Pc) - (4) 

are dipicted in Fig. 3. They are directly obtained from the data points shown in Fig. 2. 

In the second part of this paper, we propose a simple model. We show how the power- 
law behavior can be understood, and how A,„,„(p.) and A„ in Eqs. (1) and (21 can be 
estimated. The model is based on geometry and self-organized criticality (SOC) [7|, |8|. As 
we shall see, both geometry and SOC contribute powers in pt to the distributions shown in 
Eqs. (1) and (2). The relevant facts and arguments are given below: 

(A) Geometry: Let us first recall the well-known observation made by Rutherford on 
large-momentum-transfer scattering, and a less-known observation made by Williams Q| 
in which the following has been pointed out: Ordinary space-time concepts are useful for 
the semiclassical description of high-energy scattering processes, provided that the de Broglie 
wavelength of the projectile is short compared to the linear dimension of the scattering field, 
and provided that the corresponding momentum transfer which determines the deflection is 
not smaller than the disturbance allowed by the uncertainty principle. Through a simple 
realistic estimation, we can, and we did, convince ourselves that all these conditions are 
indeed fulfilled for Au+Au and p+p collisions at ^/JNN > 130 GeV and > 2 GeV/c. 
Furthermore we note that the corresponding phase-space factors can be estimated by making 



use of the uncertainty principle in accordance with Refs. 



Q and 1^ 
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FIG. 2: The power-indices, XAuAuiPc) and Xpp, evaluated by measuring the slopes in Fig. 1 are 
plotted as function of pc- 

(B) SOC: It is well-known that approximately 50% of the kinetic energy of every fast 
moving hadron is carried by its gluonic content and that the characteristic properties of the 
gluons, in particular, the direct gluon-gluon coupling prescribed by the QCD Lagrangian, 
the confinement, and the nonconservation of gluon numbers, can and should be considered 
as more than a hint that systems of interacting soft gluons are open dynamical complex 
systems which^are far from thermal and/or chemical equilibrium. Taken together with the 
made by Bak, Tang, and Wiesenfeld (BTW), these facts strongly sug gest 



observations 



the existence of SOC and thus the existence of BTW avalanches in gluonic systems 
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According to SOC, a small part of such BTW avalanches manifests themselves in the 
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exper- 



form of color-singlet gluon clusters Cq, and that they can be readily examined 
imentally in inelastic diffractive scattering processes Il3|. This is because the interactions 
between the struck Cq and any other color singlets are of Van der Waal's type which are much 
weaker than color forces at distances of the order of hadron radius. In fact, in order to check 
the existence and the properties of the Cq's, a systematic data analysis has been performed 
[l^ . the result of which shows that the size distribution Ds{S), and the lifetime distribution 
DxiT) of such Cq's indeed exhibit power-law behavior Ds{S) oc S~^, Dt(T) oc T^'^, where fi 
and u are positive real constants. Such characteristic features are known as "the fingerprints 
of SOC" P,!^. These fingerprints imply in inelastic diffractive scattering, the size S of the 
struck Cq is proportional to the directly measurable quantity xp, which is the energ y fr action 
carried by "the exchanged colorless object" in such processes, the existing data show 
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G. 3: The power-index X{pc) defined in Eq. (4) plotted as function of pc- Data are from Refs. 
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12| |. we were able to check — and only 



Ds{xp) oc X where /i = 1.95 ±0.12 [ll|,_ _ 
By considering inelastic diffractive scattering jlll . 
able to check — the existence and properties of the color-singlet gluon clusters. But, due to 
the observed SU(3) color symmetry, most of such gluon clusters are expected to be color 
multiplets which will hereafter be denoted by c*'s. Furthermore, in accordance with the 
experimentally confirmed characteristic features of the BTW theory, the SOC fingerprints in 
gluon systems should not depend on the dynamical details of their interactions, in particular, 
not on the details about the exchanged quantum numbers in their formation processes. This 
implies that Ds{S) and Dx{T) of the c*'s are expected to have jl^ not only the same power- 
law behavior but also the same power as that of their color-singlet counterparts observed in 
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inelastic diffractive scattering processes 

The fact j3| that quarks can be knocked out of protons by projectiles whenever large- 
momentum-transfer between projectiles and targets take place, has prompt us to propose 
jl^ that c*'s can also be "knocked out" of the mother proton by a projectile provided that 
the corresponding transfer of momenta is large enough where the knocked-out c*'s may 
have "color-lines" connected to the remnant of the proton. What we show now is that the 
observed power-law behavior in Eq. (2) can be quantitatively described by the product of 
the probability distribution(s) of the knocked-out c*'s and the phase-space factors associated 
with the knock-out processes. 

Recall that processes of inclusive high-pr jet-production, p+p^jet+X, at high energies 
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are dominated by two-jet events; and that in a SOC-based model such jets are produced 
in two-step-processes. In Step 1: A quark q {q^ or q^ or q^) in one of the colhding nucleon 
colhdes with a quark q {q^ or or g^) in the other nucleon where an amount of pr is 
transferred in the plane in which the two nucleons in form of thin contracted objects meet 
each other, and in which large-pr quark-quark scattering takes place. In Step 2: Since the 
two scattered g's and/or g's are in general space-like (because of the large pr), the easiest 
way for them to escape the confinement region is each "catches" a suitable time-like (in 
order to provide the high-p-j- g or g with sufficient energy) anticolor BTW-avalanches, c*'s, 
which in accordance with the SOC-picture exist in abundance in their neighborhood. This 
is how a color-singlet jet is created. A scattered g or g can also combine with a colored c* 
to form a jet or a fan which is connected with other colored object (s) through color-lines. 
This is how color multiplet jets (or fans) can be produced. Hence, in the proposed picture 
the invariant cross section Ed^a/dp^ is expected to have the following factors. 

(i) A phase-space factor that describes the chance for the two quarks (g or g, ■ ■ ■ ) which 
initiate Step 1 to come so close to each other in space that they can exchange a large pt 
(^ Et). This phase-space factor can be estimated by making use of the uncertainty principle 
and the two observations mentioned in (A) above. By choosing the z axis as the collision 
axis, Pt will be in (or near) the xy plane. The chance for two constituents (g or g, ■ ■ ■ ) 
moving parallel to the z axis to come so close to each other in the xy plane such that an 
amount pt can be transferred is approximately proportional to the size of the corresponding 
phase space AxAy ~ (Apx)^^{A.py)^^ ~ Pt^Pt^ = Pt^- 




(ii) Since each jet is associated with a to be knocked-out gluon cluster c*, which has 
energy comparable with Ej^, Ed^a/dp^ is expected to be proportional to the square of the 
probability Ds to find such a c*. The size S of BTW avalanches is directly proportional to Xp, 
thu^^oportion^ to the energy Er tt catHe. Futthetn,„.e, .tree Er « Pr for high-et>ergy 

jets [14], we have: 



This means, we expect to see a factor p^ ^ in the invariant cross section Ed^a/dp^. 

(iii) Having in mind that the scattered quarks are in general space-like, two-step-processes 
are expected to take place only when there are suitable c*'s in the surroundings immediately 
after the first step. The probability of having sufficient c*'s around, wherever and whenever 
two constituents meet during the p-|-p collision, is guaranteed when the c.m.s of one proton 





(5) 
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v=o{Pt\p + p) oc pj, ■ (6) 



meets that of the other. Hence, phase-space considerations w.r.t time requires a factor 
{Atf ~ E^^ ~ 

Hence, for p+p coUisions, we expect to see Ed?N/dp^ oc Ed?a/d?p oc p^^~'^~'^^ . By taking 
the lower hmit of /i, we obtain: 

1 £n 

2tcpt dpTdrj 

which is in reasonable agreement with the data 

Next, we focus our attention to the emj»znca/ result described by Eqs. (3) and (4) together 
with Fig. 3. Note that according to Eq. (3), the quotient Q{pt,Pc) stands for the chance 
to find a large-pr charged hadron in Au+Au collisions within a given pc range; and this 
chance is measured in "units" of the chance in finding a similar large-pr hadron in p+p 
collisions. We now take a closer look at the straight lines, on which the data-points lie in the 
log- log plots of Fig. 1, and note the fact that the slopes depend only on pc — independent of 
y/SNN- This has to be considered as a strong hint at the possible distinguished role played 
by geometry in describing/understanding such collision processes. Hence, it is useful not 
only to recall the facts and the arguments mentioned in (A) and those in (i) above, but also 

nn 

to recall that the word "centrality" pc is in fact very much involved: Experimentally jjj, |3|| 
it is determined by measuring the multiplicities of produced hadrons; but, since the notion 
of departure from the center- on- center case" is geometrical, it is expected to be describable 
in terms of geometry, in particular in terms of impact parameters, 6's. These facts and 
arguments have led us to propose the following picture. High-p^ jet-production processes 
in relativistic heavy-ion (AA) collisions can be viewed as an ensemble of corresponding jet- 
production processes in binary nucleon-nucleon (NN) collisions. The observed effects depend 
significantly on collision-geometry. 

Since every AA event corresponds to a 6-parameter (note that the reverse is not true), 
an ensemble of collision events corresponds to a set of 6-parameters. By choosing the 2;-axis 
as the collision axis where the centers of the two colliding nuclei are located at (-b/2,0) 
and (b/2,0), on the x-axis of the xy-plane in every event, we obtain point sets of impact 
parameters. It is on such point sets, the geometrical support, we study the above-mentioned 
PT-distributions. Having the observations made by Rutherford .9] and Williams in mind 
[cf. (A) and (i) above], the relation Ax ~ (Ap^)~^ ~ p^^ obtained by using the uncertainty 
principle tells us the following. For every measured value pt, there is an interval Ax; and it 
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is with this precision in the corresponding spatial coordinate that the probabihty Q{pt,Pc) 
[precisely speaking Q{Ax,Pc) on its geometrical support] of finding high-p^ charged hadrons 
can be measured. 

In the proposed picture based on SOC and geometry, we are not (at least not yet) 
in a position to make predictions for Q{Ax,pc) or its geometrical support — not even the 
dimensions of such object! But fortunately, we know how to measure them! Thank the 
master-mathematicians: K. Weierstrass, G. Cantor, H.von Koch, F. Hausdorff, Levy 
and B. Mandelbrot we learned how to use the box-counting technique. Due to the 
facts and the arguments given in (A) and (i) above, we know that the length of the boxes 
in our case are of the order Ax ~ (Apx)^^ ~ which implies: Ax becomes smaller and 
smaller for larger and larger pt- Hence the observed pr-scaling tells us that the result of this 
box-counting is nothing else but the result summarized in Eq. (3) which can also be written 
as Q{Ax,Pc) oc (Ax)'^*^^'^-'. Since this observation is independent of the positions of the boxes 
in each pc-bin, by normalizing the probabilities Q{Ax,pc), the number of boxes, A^(Ax,pc), 
needed to cover the produced hadrons distributed on the geometrical support is proportional 
to the inverse of Q{Ax,pc). That is: N{Ax,pc) oc (Ax)^'*'^^'^''. Hence, in the limit of large 
Pt, thus small Ax, X{pc) is the corresponding fractal dimension of the geometrical support 
which consists of set of impact parameter points within each given Pc-hin. This means, in 
the proposed model, we expect to see that the inclusive invariant p^-distributions for Pt > 2 
GeV/c in any kind of relativistic heavy- ion (AA) collisions satisfy 

^f^\v=oiPT\AA;p^) oc p-^--^(^=) (7) 
zTipT dpTdrj 

where Xnn ~ 7.66, and the A(pc)'s are the fractal dimensions of the point-set of impact 
parameters. 

It would be exciting to see further experiments with larger p^-values, at higher energies 
and with other kinds of colliding nuclei. Such experiments will not only serve the general goal 
of QGP-search, but also, in particular, be able to check whether the empirical regularities are 
indeed as general as the existing data seem to suggest, and to check whether/how concepts 
and methods of Complex Sciences in particular those borrowed from Fractal Geometry and 
SOC are helpful in understanding Relativistic Heavy-Ion Physics. 
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